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Abstract 

We compute the fluctuations of the magnetization and of the 
multi-overlaps for the dilute mean field ferromagnet, in the high 
temperature region. The rescaled magnetization tends to a centered 
Gaussian variable with variance diverging at the critical line. The 
rescaled multi-overlaps also tend to centered independent Gaussian 
variables, but their covariance remain finite at the critical line. 

Key words and phrases: Mean field, dilute ferromagnet, interpolation tech- 
niques, convexity methods 

1 Introduction 

The mean field dilute ferromagnet is a very interesting model. It has direct 
connections with random graph theory, as its zero temperature behavior 
reveals properties of the underlying random graph [5l [2] ■ The dilute struc- 
ture makes the model an intermediate one between the finite dimensional 
Ising model and the fully connected Curie- Weiss model. The ferromagnetic 
interactions simplify the search for basic results (such as the computation 
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of the free energy) as compared to the spin glasses, but its physical be- 
havior is still quite rich. A complete control of the model is probably as 
difficult as in the case of spin glasses and intriguing connections between 
ferromagnets and spin glasses recently emerged, so that a general theory 
of these classes of spin models is now a fundamental open issue . 

Despite its importance, the model has been basically neglected. The 
powerful cavity methods introduced by physicists could be easily guessed 
to provide the correct free energy at all temperatures and other relevant 
information, but no detailed study of the thermodynamics could be found in 
the physical literature. The community of mathematicians instead mainly 
focused on pure graph theory or on spin glasses, for which numerous results 
from physicists were followed by major mathematical breakthroughs 
while the dilute ferromagnet has been neglected. Few years ago the interest 
of mathematical physics, previously focused on mean field models, extended 
to dilute models, but surprisingly only disordered systems were studied. 
The dilute ferromagnet gained some attention only very recently. A first 
study, regarding the high temperature and the zero temperature behavior, 
was based on interpolations and convexity methods [5], in an attempt to 
complete the picture initiated with spin glasses [1] , then extended to dilute 
spin glasses [3] and to the fully connected ferromagnet as a test-case [71 
[4] within a general "structural" approach. Surprisingly, the dilute mean 
field ferromagnet has not been fully framed so far in the context of these 
"structural" methods. While the low temperature physical behavior of 
the model has not been well understood as yet, mainly due to the lack 
of physical literature about it, an important rigorous confirmation of the 
cavity ansatz for the free energy at any temperature was recently published 

Here we compute the fluctuations of the rescaled magnetization (and 
of the multi-overlaps) in the high temperature region, showing that such 
fluctuations diverge at the critical line (this holds for the magnetization 
only, as the fluctuations of the rescaled multi-overlaps remain finite) . 
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2 The model, preliminary facts, previous re- 
sults 



In this section we introduce the model and the notations, provide some 
useful formulas which are at the basis of almost all the calculations needed 
in this article, and report previous results. More details can be found in 



2.1 Definitions 

Given a set of N points, the model is defined through configurations a : 
i ~* ±1, i = l,...,Af of Ising spins. By {i^, ji,, ki,,l^}, G N, we will 
denote families independent random variables all uniformly distributed on 
1, . . . , iV. The Hamiltonian is the random function defined by 



where K is a, Poisson random variable of mean aN, for some given a G 
which is called connectivity. The expectation with respect to the random 
choice of the spins and with respect to the Poisson random variable is de- 
noted by E, and it is called quenched expectation. Given a non-negative 
real number /3, whose physical interpretation is the inverse of the temper- 
ature, the function 



is called pressure, and — Ajv(q:,/9)//3 is the free energy. Given the simple 
relation between the two, we will indifferently use either one or the other. 
The sum 



K 




(1) 



Ajv(a,/3) = -Eln^exp(-/3i7w(a)) 
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is the (random) partition function, and the Boltzmann-Gibbs expectation 
of an observable O : a C'(ct) e M is 



Zn 

a 

When it is not confusing, we will omit the dependence of on or on the 
Poisson random variable appearing in the Hamiltonian. When we omit the 
index N in the pressure we mean to have taken the thermodynamic limit: 
A{a, (3) — limAT^oo (a, P)- The main physical quantity in this model is 
the magnetization of a configuration 

N 



i=l 

often simply denoted by m. A further notation is (•) = E51(-). 

Lastly, notice that according to our notations random variables with a 
subindex zero are independent of those appearing in the weights O consists 
of, e.g. in the expression fl{aigaj„). 

2.1.1 Useful formulas 

Given a generic Poisson variable K with mean C, whose expectation is 
denoted by E, it is easy to verify that 

Ea^^ = e-'^^^'"^ (3) 

for any positive real number a. 

Another simple formula, E[(7(i^)] = ^[g{K + 1) — g{K)] for a given 
function g : N ^ R, is used to obtain 

dA^ia, 13) ^ ^^^^ /3 + E ln[l + n{cri„ ) tanh /3] . (4) 
da 

2.2 The infinite connectivity fimit 

Let us mention in what sense our model is a diluted version of the Curie- 
Weiss one, and how the latter is recovered in a suitable infinite connectivity 
limit. More details can be found in [5]. 
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Recall that the Hamiltonian of the Curie- Weiss (CW) model is 



and the associated pressure will be denoted by A''^^\(3). It is well known 
that if we let a — > oo, /3 ^ with 2atanh/3 = /?' kept constant, the 
pressure of our model tends to the one of the CW model, i.e. Ajv(q!,/3) 
A^^^\/3') uniformly in the size of the system. 

A simple proof can be obtained through interpolation 5 , considering, 
for t e [0,1], 

Ki 

N 



1 r 1 

AN{t) = -Eln^exp a,„(Tj„ + (1 - i)/?'^ W 



(5) 



i/=i 

Ki being a Poisson random variable with mean taN. 

2.3 Symmetry breaking and dilution 

Let us define 

Asia, 13) In 2 + a cosh /3 . 
We know from [5 that if 2a tanh /3 < 1 then 

A{a,l3)= lim AAr(a,/3) =ln2 + In cosh/3 = is(a,/?) • (6) 

Interpolation ([5]) allows for another simple but fundamental observation: 
applying the fundamental theorem of calculus to ^ ones observes that 
the difference between the limiting pressure A{a, (3) and its symmetric ap- 
proximation As{a,(3) decreases with the dilution. This means that such a 
difference is larger for the CW model than for our diluted version: 

AN{a,l3) - As{a,(3) = AN{a,l3) - ln2- In cosh /3 

< A^S^^ i(3') ^ In 2 ^ A^^^^ (/?') - if-) i(3') , (7) 

where we are always assuming the constraint 2a tanh /3 = f3'. This is an- 
other way the crucial difference between ferromagnetic and glassy models 
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reveals itself. In fact, in the case of symmetric distribution of the interac- 
tions the odd terms disappear in the expansion of the pressure 3 , and the 
first correction to the zeroth term (the one without spin contribution, cor- 
responding to the high temperature symmetric pressure) is quadratic and 
negative [3]. This makes the true pressure smaller than its symmetric ex- 
pression and the dependence on the order parameters (the multi-overlaps) 
convex. Moreover, the dilution increases the difference between the pres- 
sure and its symmetric expression |10j , hence when the symmetry is broken 
in the fully connected model it is broken in the dilute model too (we are 
assuming the proper temperature rescaling). In the case of ferromagnetic 
interactions, in the expansion of ([4]) the odd terms contribute [5 , the mag- 
netization being the first, and with positive sign; so the convexity just 
mentioned is replaced here by the concavity in ([4]) with respect to the or- 
der parameter Q{aigajg). This means that the pressure gets larger when 
the interactions begin contributing, and it is thus larger than its symmet- 
ric counterpart. Moreover, as explained by ([7]) the dilution decreases the 
difference between the pressure and its symmetric expression, hence the 
occurrence of symmetry breaking in the fully connected model does not 
imply the breaking of the symmetry in the dilute one. So the idea is to 
study the model at temperature zero, where it is simpler, and a symmetry 
breaking would imply the same transition at all temperatures (we again 
intend to keep 2atanh/3 — (3' constant). This is carried out in [5j. 

2.4 The transition 

We here summarize some results obtained in [5j, adding a few comments. 
The main observation is that our dilute model is in some sense delim- 
ited by its zero temperature limit and its infinite connectivity limit, both 
fully controlled [Fl [7] . We are always assuming to move along the lines 
2atanh/3 = fi' . We are about to see that the two bounds squeeze at the 
critical line /?' = 1. 
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As anticipated in the previous subsection, the study of the criticahty of 
our model cannot take advantage of the study of the criticahty of the fuUy 
connected one only, contrarily to the case of spin glasses. It is necessary 
control what happens when the dilution is decreased, rather that increased, 
and as a consequence the temperature is decreased as well. Consider 

^'^(a,/3) EE A7v(a,/3) - As{a,l3) = Apfia,!]) - In 2 - a In cosh /3 . 

It is proven in [5j that 

liminf A'^(a, /3) > if 2atanh/3 > 1 . 

N—>-oo 

The idea of the proof is to keep 2Q!tanh/3 — (3' constant, so that obviously 
2(tanh/3)dQ! + 2a(l -tanh^ = 0, and to study d A' [a, (3 [a)) / da which 
turns out to be non- negative increasing and convex in r2(cri(j(Tj(,) tanh/3. 
Such a quantity Q.{ai„aj„) measures the difference between the pressure and 
its high temperature symmetric expression, and its quenched expectation is 
obviously the squared magnetization. The model is fully solved at temper- 
ature zero |5], where it exhibits a transition, which is therefore reproduced 
at all temperatures along the critical line /?'(a, /3(a)) = 1. The convexity of 
the derivative oi A' means that along the critical line Q.{cri„crj„) takes values 
between the values it takes at /?'(l/2,oo) = 1 and at /3'(oo,0) — 1, corre- 
sponding to the zero-temperature case and the CW model respectively. As 
in both this limiting cases the model is fully solved and Er2(aioCrj„) — {m?) 
has the same behavior (the same critical index) , this means that the critical 
index is the same along the whole critical line (i'{a,f3{a)) = 1. 

3 Fluctuations 

In this section we will compute, in the high temperature region, the fluc- 
tuations of the rescaled magnetization /i = ^/Nm, which is shown to tend 
to a centered Gaussian with variance depending on the connectivity and 
on the temperature in such a way that it diverges as the critical line is 
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approached. The method extends to all multi-overlaps, whose variance re- 
mains finite at the critical line. The results are therefore very similar to 
those found in dilute spin glasses [lUj . 

The strategy we employed is the one developed in [T^ , which was then 
extended in [9l [TOj . The first step one has to take is the control of the 
way the magnetization goes to zero with the size of the system in the high 
temperature region, using a suitable perturbation. The result is that the 
rescaled magnetization remains finite in the thermodynamic limit. This 
is proven in the next subsection. This result is already desirable in itself, 
but our ultimate purpose is to prove that the distribution of the rescaled 
magnetization is Gaussian, and the proof relies on the first result. 

3.1 Bound for the overlaps 

Recall the name we gave the symmetric pressure: 

As (a, /3) = In 2 + a In cosh /? . 

We want to prove that the squared magnetization vanishes in the thermo- 
dynamic limit not slower than the inverse size of the system. Let us recall 
the definition of overlap qi...n among n configurations . . . , cr'-"-': 



For n = I one clearly recovers the magnetization m defined in 
Theorem 1 In the high temperature region defined by 



where n is a natural number larger than one and m is the magnetization 
defined in (0). 




2atanh/3 < 1 



the following result holds 



{qlJ = E17"(a,;„a,J < - (m^) = 0{1/N) , 
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Remember that io, Jo are independent of the inside f2, where v > 1. 

Proof. We will often omit the dependence of the various pressures on a 
and (3. 

We know from ^ that 

OcAn = daAs + E ln(l + fl{aigaj„ ) tanh /3) < In cosh /3 + (m^) tanh (3 

since clearly ln(l+a;) < x. But we also have ln(l+a;) > a; In 2 for < x < 1. 
Hence 

< (m^) ln2tanh/3 < c)„(^jv - As) < ^(m'^)2ta.nh (3 . 

It is clear then that we have to obtain an estimate for 9q,(Ajv — As). As in 
[ElinillO], we will not compare directly An and As- We will rather compare 
As to a perturbed pressure with a larger weight given to configurations with 
non-zero magnetization. To the purpose, let us define, for A > 0, 

A{X) = ^Eln^exp(-/3Ar(cr) + XNm^/2) 

a 

which is convex in A and such that A{Q) ~ A^- The idea is that if the field 
forcing the magnetization to be strictly positive is not too strong, then the 
magnetization will still vanish, provided the temperature is high enough. 
Notice that 

(9AA)(0) = i(m2)>0, 
which by convexity means A{X) > A{0) ^ Ajsf. Convexity also implies 

h{m^) = X{d^A){0) < A{X) - A{0) = A{X) - A^ • 

Therefore 

< dUAN - As) < i^^[A(A) - An] . 
Since 2a tanh /3 = < 1, we can choose Aq such that 

/3' = 2atanh/3 < Aq < 1 . 
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A(A(q))U=o = ^ln^exp(AoiVmV2) • 



Now let us estimate A{X) when A is chosen to depend on a according to 
A = Ao — 2atanh/J > 0, which also means A < 1. A simple calculation 
gives 

daA{X{a)) =lncosh/3 + Eln[l + OA(o-ioCTjo)tanh/3] + l{m'^)\^ 

2 da 

< In cosh /3 + (m^)Atanh/3 — {m^)xta,nh(3 = In cosh /3 , 

where the index A reminds that the Boltzmann-Gibbs measure is now de- 
fined in terms of the new weights with the field A, due to the derivative. 
Integrating back against a one then obtains 

A{X{a)) < A{X{a))\a=o +alncosh/J , 

with 

N' 

a 

Notice that this is the pressure of the CW model at inverse temperature 
Ao, which incidentally is larger than /3' but smaller than one, which is the 
critical point of the CW model. One can thus take advantage of the known 
properties of the finite size corrections to the pressure of this model, or 
simply just estimate the right hand side directly using standard techniques. 
Let us check that the finite size correction to the well known limiting value 
In 2 is of order 1/N: 

^(A(a))|„=o =ln2 + 0(l/7V) . 
Introducing the centered unit Gaussian variable J, we have 
In ^ exp(Ao A^mV2) = ;^ In Ej exp( J y^^™) 

a a 

1 ^ 1 / ^ ^ 

= InEj [] ^ exp( J VAo/A^tTi) = - InEj ( 2cosh(Jv%7^) 



1 C A 

= ln2 + — / exp[7Vlncosh(^V'V^) - I2\ 
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where Ej is clearly the expectation with respect to the random variable J. 
Now we perform the substitution y = z/^fN and use the simple inequality 
2 In cosh(y-yAo) < 2/^-^0 for Aq < 1 to get 



lln^exp(AoiVmV2)<ln2+l j exp[-iV(l - Ao)2/V2] 

= ln2+l / exp[-0V2]=ln2+l('lln— V + l 
so that 

-^ln^exp(AoAfTOV2) < In2 + 0(1/A/') . 

Notice the crucial role played by the choice Aq < 1 here. So we have now 
learnt that 



A{\{a)) < ln2 + alncosh/3 + 0(l/iV) =As + 0{l/N) 

and in the end 

„ , , ~ , 2tanh/3 - , , 2tanh/?„, , , 
< da{AN - As) < — Y^i^s - An) + — y^O{l/N) , 

where recall that As — A^ < 0. Hence it must be 

da,{AN - As) = 0{1/N) , 

which proves the theorem since, as we notice at the beginning of the proof, 
{m^) ln2 tanh l3 <da{AN - As). □ 

Let us proceed a bit further with a few observations summarizing our 
findings about the relative size of An, As, A{X). As An and share the 
same value at a = the fundamental theorem of calculus ensures that 

An-As = 0{1/N) 

and we also found that 

A{X)=As + 0{1/N) = An + 0{1/N), 0<A<1, 
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which means that, as As does not depend on A, the function A{X) has 
a very smaU variation, and by convexity this means that its derivative 
(at zero, where the derivative is smahest and gives {m?)) is also as small, 
namely of order 1/iV. 

3.2 Fluctuations 

In this section we find the probability distribution of the rescaled magne- 
tization /i — ^Jm in the high temperature region. The result is that such 
a probability distribution is a centered Gaussian with variance diverging 
as the critical line is approached. A key role in the proof is played by the 
finiteness of the rescaled magnetization in the thermodynamic limit proven 
in the previous subsection. As the results of the previous section hold for 
all multi-overlaps as well, the results we are about to prove for the magne- 
tization also extend to multi-overlaps, although their covariance does not 
diverge at the critical line. Not surprisingly, the strategy is based on the 
calculation of the generating function, following again, like in the previous 
subsection, the technique developed in [12^ and successfully employed in 

MM- 



Theorem 2 Let /i — \/Nm, where m is the magnetization defined in (0). 
Then, as N ~f oo, the variable /i tends in distribution to a centered gaussian 
variable with variance 



1 - 2q; tanh /3 

Proof. We will employ the standard method, relying on the characteristic 
generating function 

<p{u) = (exp(iu^)) , 
which we will show to be such that 



2(1 - 2a tanh /3) 
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in the thermodynamic hmit. 
It is obvious that 



du4>{u) = i{fj,exp{iufj.)} = iy N{aN eiq){iufj,)} 

by symmetry with respect to permutations of spins. 

The strategy of the proof consists in estimating the effect of the removal 
of the last spin, in the spirit of the cavity method [12l [21 UHl E] • To the 
purpose, let introduce the notations 

u- = 1/iV , ^_ = ^ cr,/\/]V^ , a_ = a(l - 1/(7V- 1)) . 

i=l 

It is not difficult to check that 

{aN enpi^iu^)) — {a n e:iq){iua n / + m_/i_)) 

and as a consequence 

du4>{u) = —u(j){u) + iVN (ctat exp(iu_/i_)) 

up to a vanishing term. 

We may now assume, with an error of order 1/N, that none of the 
terms, labelled by v and summed up, in the Hamiltonian ([1]) is a^CN, so 
to have 

,n_[exp(iu_^_)iX]<Tjv '^^ exp(/3crjv X)!^ J] 



where fi_ is the Boltzmann-Gibbs measure associated with the system with 
iV — 1 spins, at connectivity a_, the variables li, are distributed uniformly 
over {1,...,A'^ — 1}, and k is a Poisson random variable of mean 2a (see 
for instance [3] for detailed calculations of this sort). All these quenched 
random variables are independent of those appearing in the weights of ri_. 
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At this point we proceed following the idea of [12] and define 



A 



il- exp(TO_/i_) ^ ^ (Tjv exp(/3gAf ^ ^g;^)/2 



B 



g-jv exp(/3o-Af CT; J/2 



cosh'' /3 



which will be used in the following trivial identity 



B 



A 



^A AB 
B~^ 




B 



B 



) 



From the identity above is should be clear that the key idea is to simplify 
the denominator B, dealing with B instead. So we have three terms to 
compute. We want to show that the first two give the same result, and the 
third is negligible in the thermodynamic limit. 

Let us define as the expectation with respect to k and E/ = Jl!^=i ^i-v 
is the expectation with respect to the random sites appearing explicitly in 
A and B only; the other quenched variables implicitly included in are 
excluded by these expectations. 

Let us start from A, whose core is 




E„2 cosh" l3^aN Yi^iA^ + <^Ncn,, tanh (3) 




Now B clearly cancels out in the fraction A/B we are computing. At this 
point the formula in ([3]) is employed and lets us obtain 




ri_ [exp(m_/i_ )sinh(/i_2Q:tanh/3/\/iV- 1)] . 
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Now this is where the result of the previous section 

sup(/x^) < oo 

JV 

becomes crucial, as it ensures that in the thermodynamic limit the sinh 
can be replaced by its first order approximation, i.e. its argument, so that 
when — > oo we have 

iVNE^ = 2Q:tanh/3(iyuexp(iu/(x)) = (2atanh/3)9„(exp(m/i)) . 
B 

Let us proceed with the slightly more involved term AB/B, focussing 
on the numerator first 



E^EiAB = 



E^B'^fl. 



1 " 



where we used a' as a second label for configurations, since in the product 
AB there are two summations over the spin configurations. Again, B^ is 
removed by the denominator and therefore 

E-g^ = (exp(iM_/i_)- 2^ fjjvx 

Ek JJ^ (1 + (TNcri^ tanh /3 + cr'j^o-'i^ tanh /3 + (^N<yi^(T'N'^i„ tanh^ (3) 



(exp(m_/x_)2 ^ cTjvX 



Ere(l + (jNtn- tanh (3 + a'j^m'_ tanh /3 + gnc^'mI- tanh^ /S)**) 



where m_ = fi-/^/N — 1 is the magnetization of the first N —1 spins and 
g_ is the overlap between two configurations of the first A'' — 1 spins. Using 
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again ^ we get 

exp[2a(CTArm_ tanh /3 + a']^m'_ tanh /3 + aj^a'j^q^ tanh'^ /3)]) 



= (exp(w_//_)- exp(2acr^TO'_ tanh/3)x 

- (TAr exp[2Q!(crAr(TO_ tanh /3 + tanh^ 

= (exp(iu_/j,_)- exp(2a(T^m'_ tanh/3)x 

sinh[2Q;(TO_ tanh/3 + cr'^Q- tanh^ /?)]) 

and again keeping only the terms surviving in the thermodynamic Hmit we 
can write 

E— : — = (exp(m_/x_)— ^ exp(2a(T^m'_ tanh/3)TO_)2Q!tanh/3 



(exp(m_/i_)— exp(2acr^TO'_ tanh /3)cr^(7_)2Q! tanh^ /? 



= (exp(m_/^_) cosh(2am'_ tanh /3)to_)2q! tanh/3 

+ (exp(m_/x_) sinh(2am'_ tanh/3)q_)2Q!tanh^ f3 

which in the Umit reduces to 
AB 

iVNE,^^ = iVTl {exp{iu^)m)2a tanh 13+ 



i^fN (exp(m_/i_)m'(7_)4Q!^ tanh"^ j3 

i(exp(m/i)/x)2a tanh/3 = (2atanh/3)(3„(/)(u) 
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as cosh(2am' tanh /3) — > 1 and VN{m'q-) — > 0. 
The last term we have to consider is B/B: 



B 1 1 

=Ef]_[^ Vexp(/3crAr Vcr,J] = Ef]_ [- V E^(l + CTArm- tanh/3) 
B IB 2 

= Eil_[i^exp(c7jv2atanh/3)] = Ef2_ cosh(rn_2atanh/3) 1 . 



'2 

<TN 

Taking into account the orders of magnitude as before 

2 



and 



Vn\E^{1 - < ^/]Vexp(2/3a)E (^1 - |) 



since S > 1 by Jensen inequality and 

\A\ < n\eyip{/3aN^(Ti^)\ < n(exp(/3K)) < exp(/3K) . 

Thus this last third term does not contribute in the limit and the final 
result is 

(1 — 2ata.nh/3)du(pu = —u(f){u) 

which completes the proof recalling that 0(0) = 1. □ 

A straightforward generalization of the previous theorem, which re- 
quires conceptually similar but longer calculations, provides the next 

Theorem 3 Let /Ui...„ = VNqi...n, where is the overlaps between n 
configurations. Then 



1 — 2a tanh" /3 
(y^i...„^fc...fe+„) = if k^O , 

{lJ,i...nHi-m) =0 if nj^m , 
in the thermodynamic limit. 
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4 Outlook 



The computation of the free energy at high temperature and at zero tem- 
perature has been followed by progresses in two directions. The complete 
characterization (in the present work) of the statistical properties of the 
model in the same high temperature region on the one hand, the rigorous 
computation of the free energy at any temperature on the other hand [2] • 
The next important step is then a good understanding of the low temper- 
ature phase, where there are unanswered questions and interesting hints 

El- 
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